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This paper starts with the investigation of the behaviour of a set of two subsystems which are able 
to exchange some internal quantity according to a given flux function. It is found that this sytem 
exhibit a bifurcation when the flux passes through a maximum and that its kind (super-critical/sub- 
critical) depends on the dissymmetry of the flux function near the maximum. It is also found a new 
kind of bifurcation when the flux function is symmetric: we call it hypercritical bifurcation because 
it generates much stronger fluctuations than the super-critical one. The effect of a white noise is 
then investigated. We show that an experimental set-up, leading to the Maxwell demon in granular 
gas, displays all these kinds of bifurcation, just by changing the parameters of excitation. It means 
that this system is much less simple as it was thought. 

I. INTRODUCTION 

The theory of dynamical systems is a general field of research which has many apphcations in different areas, 
inchiding fluid dynamics 0, population dynamics or chemical reaction j3| among others. The idea of this field 
is to connect physical problems which are completely different and which pertain in different areas, just because 
they arc driven by a similar set of equations. Due to this, a solution obtained in one case can be applied to the 
other domains, just by settling the analogy. But to be correct, the method requires to settle the correct analogy, 
i.e. the correct change of variables and boundary conditions, which needs in turn some rigour. One of the key issue 
in these problems is the determination of the stationary states, their stability (attractors) , the dimension of the 
stable subspacc ... Also one can ask how the solutions evolve when applying different constrains on the system. 
This means to study the variation of the set of attractors as a function of external conditions to determine some 
possible bifurcation. The way to proceed is well known and it is possible to characterise the transition between 
different stationary states thanks to the bifurcation theory |H IS - This latter allows to describe the dependence 
of a stationary state with respect to some experimental parameters for instance. It has then been pointed out that 
a single stationary state may evolve to two or more stationary states as a control parameter changes. This has led 
to study the well-known (super-)critical and sub-critical bifurcations for instance. However, it is also possible that 
a stationary state evolves and splits into a continuous set of stationary states but this has not been studied yet. 
Nevertheless, as we demonstrate in this paper with a specific and simple example, this latter kind of bifurcation. 
We will even show experimentally a simple example where all these different types of bifurcations can occur, i.e. 
super-critical and sub-critical, plus the hypercritical one. In this case, this last one will occur when the system jumps 
from one kind of bifurcation to the other kind. We will firstly present the theoretical approach. Then we will present 
the classical experiment of the "Maxwell demon in granular gas" . And we will prove that this experiment can be 
used to study all these three type of bifurcation. In turn, it demonstrates that this experiment is much less simple 
than it was thought, and that the published explanations are not complete. 

We will start the paper with the theoretical analysis. To simplify the concepts, we consider a system made of two 
point-like subsystems exchanging a part of some physical quantity. This quantity is an internal property of extensive 
nature. And we investigate its dynamics. We will focus on a case where the exchange is driven by the flux functions 
of each subsystems. Each flux function will be assumed to depend only on the content of each subsystem (i.e. it does 
not depend on the other subsystem). 

When both subsystems are submitted to the same external conditions, one expects that both flux functions are 
identical. A steady state corresponds to an equilibrium state if it is stable. When the flux function has a maximum, 
multiple equilibrium states can exist. So, a bifurcation occurs when the control parameter allows the flux to evolve 
and pass through the maximum. The kind of bifurcation depends drastically on the shape of the flux curve. We will 
see that a new kind of bifurcation, called hypercritical, can be found with some special constraints on the shape of the 
flux function. It corresponds to the case when a unique attractor degenerates into a subspace of dimension larger 
than 0. This generates an indeterminacy at some precise value of the constrains. Apart from this precise constrain, 
cither a single attractor or a set of two attractors exist. The case which will be studied here is the case = 1 but 
larger dimension of subset could be envisaged (d^ = 2, 3) or fractal subsets too d-A > 0. Most of the analysis will 
use a continuum approach. However, at the end of the theoretical part the effect of noise will be introduced; this 
one can be due to the "granular" nature of the extensive quantity. When this noise is taken into account, the kind 
of bifurcation is particularly ticklish to display and the set of the two subsystems undergoes a diffusion-like process. 
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which may be solved using a Langevin's formalism. This diffusion process may broaden the "point-like" attractor 
when the attractor is unique, or when the attraction of the attractors is large compared to the noise agitation when 
there are more than one attractor. In this case the steady state expands over a given volume of the phase space. 
But the noise can also generate repeated jumps between the set of attractors when attractors are discrete in this last 
case, or even provoke a diffusion among a continuous or discrete set of possible attractors. All these configurations 
will be considered. 

In the last part, the paper will revisit the problem of the Maxwell's demon in granular matter, which can be stated 
as follows: consider two halve-containers containing few grains each and connected by a lateral hole. The system is 
vibrated and the grains are agitated so that they can pass from one container to the other one alternately. In practise 
one should expect an equi-repartition. However this one is ensured only when the vibration excitation is large enough, 
but when the excitation becomes two small and the dissipation is large enough, one observes that a container is more 
filled than the other. This occurs because the larger dissipation in the more filled container reduces the agitation of 
the grains and the probability of these ones to escape, while the smaller dissipation in the second container increases 
the speed of the grains in this box and favoured the transfer to the other box. The present experiment will use a 
single vibrated container containing some grains, and in which a hole has been drilled on a side to allow grains to 
escape. This will allow to measure the flux function J of grains going out from the box as a function of the number 
N of grains in the container. We will see that the measured flux functions let predict the existence of an hypercritical 
bifurcation for some range of parameters. As this result was not found by previous experiments, numerical simulations 
and theoretical approaches, it proves that the problem was not properly settled and that the correct equations are 
still not given. This proves that much work remains to be performed on granular materials and that much care has 
to be taken to settle correct analogies. 

The paper is built as follow: Section ^1 presents the general issue, the equations which rule the evolution of the 
system and the equilibrium condition for which the internal quantity in each subsystem remains constant. Section 
mil presents an analysis of the bifurcation problem occurring when the internal quantity maximises the flux function. 
The hypercritical bifurcation is defined and the condition of its occurrence characterised. The effect of noise is 
also discussed. Finally, section Hvl is devoted to a convenient experiment to study the so called "Maxwell's demon 
phenomenon in granular matter" presented above. The paper is concluded in section Ivl 



Consider two subsystems, left and right (I, r), each characterised by an internal extensive quantity, xi and Xr and 

submitted to a set of external conditions, and j/r . We assume first that the set of the two subsystems is closed 
such that: 



Also the exchange between the two subsystems is ruled thanks to two flux functions, Ji and Jr associated to each 

(i) (i) 

subsystem. These flows are assumed to depend on the external parameters, j/; or ?/r , applied to each container 
and on its content xi or Xr only. Thus, this paper aims at describing the evolution of such subsystems, at finding if 
stationary or equilibrium states exist, if they are stable (equilibrium states) or not and at studying the possibility of 

exhibiting some bifurcations when external parameters, or yi^\ evolve. 

As an example, one can consider the "granular Maxwell demon" 0, where two containers connected by a slit 
are vibrated. These containers are partially filled with macroscopic particles and can exchange some particles across 
the slit. Here, the internal quantities are the number of particles in each box while the external conditions are the 
dimension of each box, the position and size of the slit, the frequencies and amplitudes of vibration for each box, 
the gravity and so on. In this peculiar case, it is known that particles equipartition is observed for intense enough 
vibration, but that it breaks at small enough vibration excitation, displaying a threshold and a bifurcation. 

Coming back to the more general case, the time evolution of the system is given by the time evolution of xi and Xr 
which writes: 



II. PRELIMINARIES 



Xr + Xi = Xtot = 2xo = cste. 



(1) 



= Jr{Xr) - Jl{xi). 



(2) 



And equation imposes: 



dt 



dxi 
It' 



(3) 
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dxi 



Both subsystems are stationary (and then in equihbrium but not always stable) if — ^ 
that their flux are equal: 

Jl{xi) = Jr{Xr). 



0, which imposes in turn 



(4) 



One has now to determine whether this steady/equilibrium state {xifl, Xr,o) is stable or not. It is stable if any small 
perturbation Sx of xi^ and Xr,o (i-e. Xi = xi^ + Sx and 6x) decreases spontaneously with time. It is 

unstable on the contrary. A first order expansion gives: 



d{xifi + Sx) 



dt 



dSx 

~dr 

dSx 

~dr 



Jr{Xrfi - Sx) - Jl{xifi + Sx), 

—Sx 



dJr , \ , dJi 

dx dx 



-aSx. 



(5) 
(6) 
(7) 



Where a 
a < 0. 



dx 



(xr.o) + ^(2'';,o)- Consequently, a steady/equilibrium state is stable when a > and unstable when 



For the sake of simplicity, we will limit hereafter the formulation to two systems with some symmetries, but 
generalisation is straightforward. So, we will assume that both systems are identical (same size, slit at the same 
position, ...) and are submitted to equivalent external conditions. It means in particular that J; and Jr depend only 
on xi and x^ respectively and on a unique set of external conditions. This imposes also that Ji{x) — Jr{x) = J{x). 
One would like now to know if a stable equilibrium state can become unstable when varying a control parameter. This 
leads to determine under which condition a bifurcation could occur and which kind of bifurcation could be obtained 
from this formalism. 



III. BIFURCATION ANALYSIS 



A. Classical analysis 

Let us assume that the two subsystems contain both xo initially. From equation Q and according to equation 
as soon as the two subsystems are identical and are submitted to the same constrains, the state defined by 
xi = Xr = xq is a steady state. According to equation Q, this state remains stable as long as the first derivative of 
J at xi/r is positive that is, as long as J increases with x at xq. This solution becomes unstable as soon as J starts 
decreasing. In this case, a new solution has to be found which satisfies xi + Xj. = 2xo = cste and J{xi) — J(a;,.) = 0. 
This is the jump from the unique solution to the set of two different solutions which will interest us during this 
section. It occurs when dJ /dx passes through 0. 

Generally, as J is an outgoing flux, one expects that J increases at small x starting from J(0) = 0. So, the 
jump/bifurcation occurs when J passes through a maximum. However, it may occurs that J starts from a non null 
value and passes through a minimum in some cases but the analysis is similar to the one made in the following when 
J increases at small x. 

Note also that we will consider in this paper only the flux function for which the first derivative could cancel at most 
one time and the case where J is infinitely derivable. 

So, turning back to the most probable case for which J increases at small enough x, i.e. ^ > at small x, it 
means that a transition occurs when the flux function reaches its maximum value J(x,„). This leads to study and 
determine the kind of bifurcation which occurs around xi^r = Xm- The parameter v which controls the distance to the 
threshold is v = xi -\- Xj. — 2xm and can be used to study the bifurcation nature. In this case, the preserved quantity 
xi + Xr writes xi +Xr = 2x,„ + v. It is also convenient to introduce the asymmetry parameter the order 

parameter. So writing J as a function of the distance to the maximum and using a Taylor expansion, one obtains: 



E 

oo 

E 

fe=i 



{Xl ~ Xm) d^J ^ 

kl dx'' 

(^Xf Xfii^ d J 

k\ dx'' 



Xm), 



(8) 
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We use the notation hereafter = ^{xm- By hypothesis, J^'^\xm) = ^{xm) < andJ^^^ (xm) = ^(x^ 



is the bifurcation threshold. 
So, using the change of variable: 
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Xl ^m; (9) 



The dynamics equations becomes at third order 

du 



0. (11) 



J'^^''{xrn) = and j(^)(a;„i) < by hypothesis. 
The solution of equations Hll() at first order is: 

u(0 ~exp-(wj(2)(x„)t). (12) 

As J^'^\xm) < 0, the solution defined by equation (|12|) is stable for v < and unstable for v > and one deduces 

the typical time t of evolution r = (vJ^'^"'{xm)) ^ ■ It is worth noting that t tends to infinity as v tends to 0. This is 
the so-called "critical slowing down" and it means that the equilibrium needs an infinite time to be reached. To get 
the correct time behaviour at w = 0, one shall expand equation Hll() to higher order in u. This leads to a power law 
relaxation instead of an exponential one. 

A stationary state corresponds to ^ = 0. Thus, a solution to equation (|ll|l with ^ = is: 

u = 0. (13) 

As already told, this solution is stable if u < and unstable if u > 0. It means that stable equilibrium 

state as long as xi + Xr < 2xm- So, when v > the solution u = is unstable. To get the behaviour, one has to take 
account for higher order term in the development. 



1. Case J'^^\x^) = 0, J^'^\x^) < 0, J^'^\xrr,) > 
Limiting to the third order, this gives two new solutions for equation 111() when J'-^-'(.t,„) > and w > 0: 



U± = ± \ —V 



1/2 



= ±VM- (14) 

One can show that these solutions are stable. As the two new solutions do not exist at w < 0, the bifurcation is 
the one described in figure ^ with a parabolic branching at v > 0. This is characteristic of a super-critical fork 
bifurcation (because the new solutions are stable). In the present case, it occurs at J^^-'(xm) = 0, with J'^^)(a;,„) < 
and J^^\x,n) > 0. 

The two solutions appear also when j(^^(a;m) > and v < but this case will be studied later. Let us simply 
mentioned that they are unstable. 

To conclude with the case {J'-^^Xm) — 0, j'^-'(xm) < 0, j'-^^Xm) > 0}, when v increases from w < to w > 0, the 
solution xi ~ Xr then breaks following a super-critical bifurcation with two symmetric and stable solutions u± shown 
on Figure n] Figure |21 shows a typical shape of J when J^^-'(xm) < and J^'^)(a;m) > 0. It is characterised by a 
left wing steeper than the right wing. This allows the x position of the middle of an horizontal secant at given J to 
increase as J decreases, which is the requirement to get the super-critical bifurcation. Indeed, in the case 2x > 2xm, 
any steady solution above Xm requires to find two solutions xi and X2, such as J{xi) = J{x2) and X1+X2 = 'ix > 2xm- 
Such a possibility exists only if the middle of the horizontal secant at a given J increases when J decreases. This is 
also what J''^^{xm) > means when J(^'(a;„0 = and J^^^x^n) < 0. 
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FIG. 1: Solutions of equations (lllll for J^^\xm) > 0. Full lines show the stable solution while dashed line shows the unstable 
ones. Arrows describe the time evolution of a perturbation 5u. Super-critical bifurcation when v becomes positive. 




FIG. 2: An example of the shape for J corresponding to a super- critical bifurcation (i.e. with j' ' (xm) > 0). The right wing 
shall be broader than the left wing. The dashed line is the location of the middle of horizontal secant. 



2. Case j'^) {xm) = 0, f (xm) < 0, j'^^ (x™) < 

When j'-^^Xm) < 0, there is no solution to equation Hll|) with ^ = for v > 0. This is due to the fact that the 
position Xmid of the middle of the secant at a given J increases with J when the right wing of the J curve is steeper 
than the left one, as it is exemplified in Figure 0] So, to get a solution one shall expand equation (jS)) at higher order 
in the vicinity of = in equation Furthermore, owing to the fact that the dynamics of the system is controlled 
by the subtraction of two flows, the w"' term cancels (as the term had disappeared in equation and one shall 
develop to fifth order. 

Using the change of variable (lO and ifTUIl ) and applying these expansions lead to the following dynamics equation: 

A possible stationary state (^ = 0) of equation p5|l is u = but it is only stable if u < in the vicinity of v = 0. It 
is unstable for t; > 0. The other possible solutions are symmetric compared to m = 0. Furthermore, equation H15() with 
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^ = accepts three solutions when v > (two stable and one unstable), five solutions in the range v S [vmin < 0; 0] 
(three stable and two unstable) and one stable solution when v < Vmin- 

A typical example of such a case is displayed on Figure 13 An example of curve J corresponding to this bifurcation 
is displayed on Figure^ It is characterised by a right wing steeper than the left wing, as explained previously, and 
by the position x of the middle of a horizontal secant that increases when J increases beyond a value depending on 
its shape. 

So, when v < Vmin, w = is the only (stable) solution. Then, when v G [vmin < 0;0], the solution u = remains 
stable but there are also two other stable solutions, both are separated from the u = solution by an unstable 
solution. This is typical of multistable state, which is known to lead to hysteresis. When z; > 0, u = is now an 
unstable solution and the two symmetric solutions u± are the two stable solutions. Hysteresis occurs as follows: 
we start with u ~ sd v < Vmin- So, increasing slowly v from v < Vmim let the solution u = unchanged, then 
increasing v above v ~ forces the jump of m = to u+ or u- when v becomes positive. Now, decreasing slowly v 
from above to below v < Vmin the solution u+ or u- evolves slowly till v reaches v < Vmin where u jumps to u = 0. 
This hysteretic behaviour is typical of a sub-critical bifurcation. 
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FIG. 3: Solutions of equations IIH for j'^^^Xm) < 0. Full lines show the stable solutions while dashed lines show the unstable 
ones. Arrows indicate the spontaneous direction of motion of a perturbation 5u. This is typical of a sub-critical bifurcation. 



J 




FIG. 4: Typical shape of J leading to a sub-critical bifurcation (with j'^^^Xm) < 0). The dashed line corresponds to the 
location of the middle of horizontal secant. 



For w = 0, that is xi + x^. = 2xm = Xtoti the dynamics equation writes: 
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The dynamics is controlled by J^^\xm) and J^^^{xm)- Return to equilibrium is no more an exponential decrease but 
a power law. The steady solutions of equation (|16|l writes at first order: 



w = 0, 
u± = ± 



1/2 



(17) 
(18) 



if, j'^^\xrr.)J^^Hxm) <0. 



Their stability depends on the third and fifth derivative of J at x^- Finally, the comparison of Figures |21 and 0] shows 
that the symmetry of J plays an important role on the bifurcation behaviour. 



3. Symmetric flux function: J{xm — x) = J{xm + x) 



A special attention has then to be paid to the case of a symmetric J, i.e. with two symmetric wings, which implies 
J^''\xm) = for all odd k. The dynamics equation writes at fourth order in such a case: 

§ = -uv[j^^\x^)+^^^j('\x^)^ 

-U^y^'\xra) + ... (19) 

Where J^^\x„,) < 0. 

Let us first neglect the higher order terms J'^''^(xm) and consider the problem near v = 0. As J(2)(.T,„) < 0, the 
solution u = is stable for w < and unstable for v > 0. 

So, at u = 0, there is a continuous set of equilibrium states, since ^ = whatever u. So, all initial states are 
equilibrium states: keeping v = 0, and starting from any definite state u — uq, uq is a steady state. Suppose now that 
the system passes spontaneously to ui due to some perturbation, in this case the new ui state also is the new steady 
state and the system will be fully driven by perturbation or uncontrolled noise that generates transition from uq to 

Ul. 

As we will see, such a bifurcation trend is important because it leads to generate extremely large fiuctuation, larger 
than a critical bifurcation. We can then call this bifurcation an hypercritical bifurcation. Its diagram is displayed on 
Figure El Figure shows typical functions J which give rise to a hypercritical bifurcation, with the vertical dotted 
line representing the location of the middles of horizontal secant. 

This bifurcation is particularly interesting from an experimental point of view because its dynamics is completely 
driven by the existing noise, and because any extensive quantity is subject to noise due to the discreetness of nature 
at some microscopic finite scale. Indeed, no matter some mechanism exists to absorb any perturbation or to make the 
perturbation grow, the system is then entirely driven by the noise when it is at this working point even if it covers 
only equilibrium states. 
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FIG. 5: Solutions of equations IIH for v = and du/dt — for all u. Full line shows the stable solution while dashed line 
shows the unstable ones. The dot-dashed line shows the set of solutions which are neither stable nor unstable. Arrows indicate 
the spontaneous direction of motion of a perturbation Su. Hypercritical bifurcation. 
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FIG. 6: Two examples of the shape of J corresponding to a hypercritical bifurcation when v — 0. The curves are symmetric 
with respect to the maximum of J, J{xm)- The dashed line is the location of the middle of horizontal secant. 



Does it means that the dynamics of the system at w = is a simple diffusion on the u axis? Not com- 
pletely, because the fluctuations that generate the motion in the u direction may have specific characteristics 
that make the system more complex. For instance, even in the case of a random process imposed by discrete 
nature of the flux, as the noise is often linked to the square root of the magnitude of the flux, one shall expect 
already that the diffusion coefficient shall depend on the value of u, because the value of J depends on u (see FigureEl. 

Since J(0) = 0, J{2xm) ~ and the fiux becomes negative when v > 0. This is why figure [S] does not exhibit 
solution for v > 0. However such solutions may appear in this hyper-critical bifurcation if some dissymmetry of J is 
introduced at large v only, i.e. above x = 2xm — a with a being small compared to Xm, so that one will find a solution 
for V > 0. This requires of course that J{x) remains positive for all x > 2xm — a and tends to as a; tends to infinity. 

4- Slightly dissymmetric flux function 

If the flux function is slightly dissymmetric in the range [0; 2xm], that is, if some of the odd derivatives at a; = Xm 
arc not null beyond order 2fc -I- 1 > 5, the bifurcation is not hypercritical stricto sensu. It means that a solution 
may be given. However this would require a large/inflnite precision especially in the vicinity of w = to maintain its 
state at the predicted value. On the contrary, let us assume some existing intrinsic fluctuations. It means that the 
system may explore different states with time, and the wider the range explored the nearer from the symmetry the 
J function. This is because the dissymmetry of the J curve is linked to the derivatives of small order {2k + 1) at the 
maximum. So, one may expect that in the vicinity of v = 0, the larger the order k for which j(2*:+i) {xm) is non null, 
and/or the smaller the J'-^'^^^-' (x„i) the wider the range of u over which the system may spread due to fluctuations. 
It is then quite important to include the effect of fluctuations in the dynamics to treat correctly the problem. This 
will be done in the next sub-section. 



B. Effect of the noise 

The previous section was devoted to the bifurcation analysis when symmetry breaks. Expected behaviour has been 
described without taking into account any noise and it has been found that the kind of bifurcation depends strongly 
on some detail of the shape of the flux function. This makes the solution sensitive to tiny details. On the other hand, 
it is experimentally impossible to completely cancel the noise. So, one has then to introduce it to describe properly 
the physics in the vicinity of the bifurcation. 

Among the questions to answer, one would like to know if the system is able to jump spontaneously from a stable 
solution to an other one just because of the existing noise. 
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This depends probably on the distance between the stable solutions and on the intensity of the noise. One can also 
ask how the system explores the space of stable configurations in the case of a continuous set of stable states, or how 
far from a stable state the system can be when it is driven by a noise of given amplitude. These questions can be 
addressed thanks to the computation of the fluctuation of x around a stable solution. But this requires also to model 
the action of the noise. 

Let us start with the last problem and consider a stable state [xi; Xr) and some perturbation 5x , i.e. [xi+5x^ Xr—Sx). 
Without perturbation, its dynamics writes: 

dxi 

--^ = J{xr)-J{xi)=0. (20) 

Let us slightly perturb this equation and use a Taylor expansion at first order: 

d{xi + Sx) 



J{xr-Sx)- J{xi+Sx), (21) 



dt 

dSx 
~dt 



dSx ^ 

^ w ~aSx. (22) 

dt 

where we have used a = j'^-'(x-() + j'^^'>{xr) > (because it is a stable state) and J{xi) ~ J{xr)- 

Let us submit this system to some random white noise of amplitude s{t), of zero mean and of variance ct^, the 
system obeys to a Langevin-like equation : 

dSx - , , 

— « -a5x + e[t). (23) 

A characteristic time is the time for the system to exchange a "particle" (or quantum) of the extensive quantity 
x: Tc = j(^xi)+J{x ) • Assuming that t is always much larger than the time for xi and x^ to exchange one particle, 
Tc ^ j{xi)+j{x ) 1 equation H23|) can be solved following the classic procedure as in the case of a Langevin equation 
and find the variance crf^ of 5x, as a function of the variance of the noise, We start from the general form of 

the solution of equation lP5|l : 

5x{t) = fa(0)e-"* + /* e-"(*-*')e(t')d<'. (24) 
Jo 

As the noise e is a zero-mean random variable, (e(t)) = 0, the expectation of 5x writes: 

{5x{t)) = fe(0)e-"*. (25) 
So, the perturbation damps, in average, from any (small) initial value with a typical time r: 

r=i. (26) 
a 

Considering now the variance of 6x around a stable position, it is defined as: 

ai = {5x{tf) ~ {5x{t))\ (27) 
Using the square of equation H24|l . one obtains: 

aj^ = f dt' r(e(f')e(i")>e~"^*"*'^e-(*'-*")di". (28) 
Jo Jo 

As £ is a white noise, one finds: 

^L = g(l-cxp(-2at)). (29) 

So, when t = 0, the variance is null because 5x{Q) is known exactly. As t increases from 0, asx starts to grow as the 
square root of time under the effect of the noise. At larger time, the standart deviation saturates and reaches <7^/\/2a.. 
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Note that, we shall only consider the case where asx is small since a = cste requires essentially &x to be small. 
So, one obtains: 

ol, - o% (30) 

if Tc ^ t <C a^^ . The width of the distribution of then spreads like \Jt as in a diffusion process until t w oT^ . As 
t 3> oT^ , the width of the distribution of 8x around a stable position is 



i. /n the vicinity of the bifurcation 

Let us assume that the stable state of equilibrium is in the vicinity of Xm ■ If xi and Xr are close to Xm such that J 
is almost constant, then J{xi) « J{xr) and a « cste <C 1, one obtains (when t <C a^^): 

al « a^t. (32) 

(5x undergoes a diffusion process and nothing prevents from the spreading of the perturbation as long as Sx is small 
enough, that is, as long as asx is smaller than the size of the region where a ^ 1. In this case, the diffusion coefficient 
writes D = cr1/2. The time needed by the perturbation to reach any small distance d from an initial position is: 

So, if two stable states are close enough, the sytem will pass from one stable state to the other with a frequency 
equal to al/2d'^. diffusion has then two effetct: the broadening of each steady state and the jump from one steady 
state to the other. 

In the case of a supercritical bifurcation, we may try to define when the two " steady" states remain separated and 
detectable as such. It occurs when the standard-deviation (due to noise) of a single steady state (without noise) is 
smaller than the distance between the two steady states xi = x^ = Xtotl'^ and Xm- As the standard deviation ugx is 
given by eauation l31l it imposes: 

2 

al,^^<{x^-xif. (34) 
With Xm — xi — —v/2 and a = 2J'^-'^^ (x//^) ~ ^vJ^'^^x^)^ one obtains: 

2vJ(^)M ^ ^' ^^^^ 
=>al < i«3j(2)(x„). (36) 

If this inequality is not satisfied, the two branches of stable solutions in figure [3 can not be observed and the 
bifurcation is not detected as long as v is too small. 

In the same way, when a subcritical bifurcation occurs, as displayed in figure^ the jump from the single stable 
solution M = to one of the two others may occur before the maximum of the flux function, i.e. before u = 0, so that 
for strong enough noise or little subcriticality, the bifurcation may look as a supercritical one. 

Finally, if one considers a hypercritical bifurcation, the flux function is symmetric with respect to its maximum, a 
is always null and J{xi) — J{xr)- All the states are equilibrium ones and, because of the noise, Xi and Xr constantly 
fluctuates between all of them. As a conclusion the nature the noise can hide the true nature of the bifurcation. 

2. Far from the bifurcation 

On the other hand, if the stable state of equilibrium corresponds to different xi and Xr such that a w cste 1, the 
variance cr|^ tends to a finite value: 

al « al/2a. (37) 
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It means that the spreading of 5x due to the noise is balanced by the damping due to a. Assuming that a is constant, 
Sx is a Gaussian random variable having the following probability density function: 



1 



a 

2^ 



(38) 
(39) 



The probability to jump from a stable solution to another one is the probability to have 5x larger than the distance d 
from the initial stable state to the frontier of its basin of attraction. As a has a damping effect, one expects that the 
probability for a jump to occur should decrease as a increases. From equation (|39|1 . one can write this probability: 



P{5x >d)^\ 



1 - erf 



\/2ad 



(40) 



where erf(.) denotes the error function and P{5x > d) decreases as a increases as expected. 

In fact, the latter result can be applied if a supercritical or a subcritical bifurcation has occured and not after a 
hypercritical bifurcation. After a hypercritical bifurcation the only result one can write is that the life expectancy, ti , 
of a stable state is proportional to the time needed by the noise to make the two boxes exchange one particle: 



1 



J{xi) + J{Xr)' 



(41) 

(42) 



So, as J{xi) + J{xr) tends to zero, the corresponding equilibrium state will be observed during a time much longer 
than the one for which J(a;/) + J(a;,.) is large. Thus, from an experimental point of view, much attention has to be paid 
to the interpretation of the experimental results. Indeed, a state may be wrongly interpreted as a non equilibrium 
one if J{xi) + J{xr) is large. 



IV. EXPERIMENTAL PART 



A well-known experiment involving a system of two communicating containers partially filled with macroscopic 
spherical particles is the so-called granular "Maxwell demon" . When such compartments are shaken together at 
same frequency and amplitude, it has been often quoted that, under some circumstances, one compartment empties 
spontaneously in the other 0. This symmetry breaking comes from the fact that the outgoing flux function is not 
monotonic as the number of particles increases: it first grows as the number of beads N increases and, as TV > N^ax, 
decreases because the dissipation becomes more important. Then, it seems to be a relevant experiment to test the 
prediction of the previous section. Note that this experiment has been studied in Q. He used the equality of the 
outgoing flux from each box as the equilibrium condition, and a continuous model taking into account that the 
effective granular temperature decreases when the density of particles increases. He concludes that the symmetry 
breaking follows a second order phase transition. However, it has been pointed out that the Egger's approach suffers 
some discrepancies with respect to experiments |lCj and that a thermodynamics analogy to describe such a system is 
questionable • An other approach has been proposed when the size of the hole connecting both compartments 
is larger than the mean free path of the particles. It is based on the pressure balance between the two boxes and 
leads to the conclusion that beyond a critical value of the number of particles, the symmetry breaks spontaneously. 
Anyway, this surprising behaviour has motivated many experiments usin g m ixtures of different particles, different 
boundaries conditions or different number of compartments for instance j3 ll^ ITEL . But the influence of the 
experimental parameters on the shape of the flux function, and then on the kind of bifurcation leading to the 
symmetry breaking, has not yet been addressed. 

In the present experiment, we study a granular Maxwell demon experiment using the outgoing flux function ,/ 
to forecast what kind of transition occurs when the symmetry breaks. The control parameter is the number of 
particles N and the analysis is made for three different frequencies. We use a single vibrated box to measure J as a 
function of the number of particles N that it contains and of the frequency of excitation. Our analysis is based on 
the same equilibrium condition than the one used in , that is on the balance of the outgoing flux between the two 
compartments. Note that because the number of balls is discrete, there is an intrinsic noise in our experiments but 
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we will see that it is nevertheless possible to identify different bifurcations. It is also important to note that in the 
case of two vibrating boxes and if N is large enough, interactions between particles coming from one compartment to 
the other could occur in the vicinity of the route between the boxes. However, this latter would occur at much larger 
densities. It would lead in this case to a flux function depending on the number of particles in both boxes which is 
not our purpose here. 



The experimental apparatus is shown on Figure |7| A single box of dimension L = 20mm, H — 30mm, I = 12mm, 
with a narrow slit on one side at height IQmm and of width Smm, is partially filled with TV w 680 steel spheres of 
diameter 1mm. 



The box is fixed on a sliding girder connected to a crankshaft and a rotating wheel. It allows to impose a linear 
oscillation to the box, A sin (27r/t) where / and A are the frequency and the amplitude of the sinusoidal oscillation 
respectively: / G [26Hz — A3Hz] and A ~ 1.25 mm. The mechanical device is mounted on a massive steel socle in 
order to avoid any perturbation. 

The frequency is imposed with an electrical motor and controlled, during each experiment, thanks to a tachometer 
to measure the fluctuaction of the frequency. The maximum amplitude of these fluctuations is smaller than 2% of 
the imposed frequency. When beads leave the vibrating compartment, a long tube bring them to a repository with 
a soft bottom such that grains do not bounce. The evolution of the mass of the content's repository is measured 
each 0.5s thanks to a weighing machine with a precision much less than 10~^g (which is the tool accuracy), because 
of parasite effect that includes mechanical vibration and ball collision with bottom. Data are then recorded with 
Lab VIEW Software on a computer and smoothened. For a given frequency /, the outgoing flux J{t) is derived from 
the temporal mass evolution. This allows to finally plot J as a function of the number of beads, N, remaining in the 
vibrating box. 



Figure IHl shows the outgoing flux, J{N), obtained for three imposed frequencies / = 2Q.lHz, 33.3i/z and 42.2Hz. 
Each plotted curve is the mean of several experiments. They all starts from J(0) = 0, increases until it reaches a 
maximum and then starts to decrease. This behaviour reminds us of the different shapes of flux function we studied 
in the previous section. 

One can already note that Figure|Hla has a left wing steeper than the right wing and that|S|b displays a right wing 
steeper than the left one. So, if the experiments are described correctly by the model of the previous section and even 
with some noise, we may be tempted to conclude that the kind of bifurcation changes as the frequency of vibration 
increases. 

Assuming a Gaussian white noise, the standard deviation (s.d.) is computed by multiplying the empirical s.d. with 
the 5%-quantile of the Student law with p degrees of liberty, where p is the number of experiments Indeed, it 
is a classical way to proceed when the number of experiments is small. As one can expect, at fixed frequency, the 
uncertainty is large when the flux is large. Note that three experiments have been made for / = 26.7Hz and nine 



A. Experimental apparatus 




weighing machine 



FIG. 7; Experimental apparatus. 



B. Experimental results 
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FIG. 8: Experimental flux as a function of the number of particles in the container shaked at frequency f. a. f — 26.7Hz, h. 
f = ii.mz, c. / = 4a..2Hz. 



experiments have been achieved for / = 33.3i?z and / = A2.2Hz. This explains why the uncertainties seem larger 
for / = 26.7Hz than for the two other frequencies used. 

One has also to mention that in these experiments and for the largest values of , the number n of layers of particles 
is n « 3. The distance between the top layer and the slit is then Z? w 7 mm. For the highest frequency, / = 42.2 Hz, 
the height reached by a particle initially at rest is ft- w 5.6 mm. As h ^ D, the density of particles in the vicinity of 
the slit might be sufficiently large such that interactions occur between the beads of the two boxes. It means that 
the tail of each curve in Figure |S1 that is for the largest number of particles, could have been modified by the high 
concentration of particles near the slit. In this case, this would involve another equilibrium condition than the one 
on flux and then an approach using fluid parameter, similar to |l2l | for instance, could be more suitable to this problem. 

Anyway, if one assumes that interactions do not occur, one can use the result of section lllll Each curve is then fitted 
to a polynomial function of degree chosen thanks to a likelihood ratio test . One then computes the value for 
which dJ/dN = and the successive derivatives of J at N^. The result is that Figures|Hla,|Slb and|Hlc correspond to 
supercritical, hypercritical and subcritical bifurcations respectively as N increases. As forecast previously, the kind 
of bifurcation depends on the frequency used during an experiment. 

Let us discuss briefly how the system converge to an equilibrium stable state and the effect of the noise. The 
convergence to a stable state depends on the sum of the derivatives of the flux function at Ni and Nr. However, as N 
is a discrete value, the flux of N is discrete and fluctuates, and the exact equilibrium state is not reachable. It results 
a standard deviation at Ni and Nr which is inversely proportional to the sum of both slopes around Ni and Nr (see 
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equation 13 1() . 

Recently, Mikkelsen et al studied a granular Maxwell demon experiment and the effect of statistical fluctuation 
on critical phenomena. Their numerical simulations show qualitatively the same kind of flux curves as ours, as the 
frequency increases. The symmetry breaking is studied with the frequency as the control parameter and with different 
total number of particles N to highlight the effect of a noise. They find that, at fixed N, the transition is supercritical 
when the frequency becomes larger than a critical threshold. However, according to our analysis, it can be either 
supercritical, subcritical or hypercritical when N becomes larger than a threshold, the type of bifurcation depending 
on the frequency of vibration. They also find that they arc faced with a supercritical bifurcation when the flux curve 
is symmetric whereas our study proves that the bifurcation is hypercritical and that the system is driven by the noise 
in this case. This shows how the type of bifurcation depends on the control parameter used. 

As a final remark, it is worth mentionning that a Maxwell demon experiment checks directly the breaking of symmetry. 
It tests then the shape of the flux function near its maximum and cannot then be considered as a proof of any modelling 
of the flux function far from this maximum. Consequently, most of the flux function assumed in the litterature cannot 
be considered as settled on experimental evidence. 

V. CONCLUSION 

In conclusion, wc have proposed a theoretical analysis of the behaviour of two subsystems able to exchange each 
other some part of their internal quantity thanks to a flux function J. We have shown that when J is not monotonic, 
a bifurcation can occur leading to different equilibrium states. The type of bifurcation is then found to depend on 
the detail of the shape of the flux function near its maximum. In particular, we have pointed out a new kind of 
bifurcation called hypercritical when the shape of the flux function is symmetric with respect to its maximum. In this 
case, a continuum of solution exists making the dynamics of the system govern by the noise. In the other cases, the 
noise, if large enough, is found to have an impact on the conclusion concerning the kind of bifurcation which appears. 
In the last part of this paper, we have presented a simple granular Maxwell demon experiment to investigate the kind 
of transition occurring when the symmetry breaks. Unlike previous experiments, we use the number of particles as 
the control parameter with different frequencies of vibration /. Our results show that different bifurcation, including 
the hypercritical one, occur depending on /. This behaviour was not pointed out before neither experimentally, nor 
numerically, not theoretically. It could then allow to investigate new kind of bifurcation. 

However the detection of the true nature of bifurcation (supercritical vs. subcritical, etc) is made difficult due to the 
little number of grains (and by the noise it induces) in the case of a Maxwell demon experiment. This is easier if 
one uses directly the flux function which gives directly the symmetry of the curve. In turn this means the use of a 
single half box. We believe that a Maxwell demon experiment with 100 balls only test the presence of a bifurcation, 
but not the precise shape of the flux function near its maximum or further in the tails. So, previous Maxwell 
demon experiments can not be considered as a test of the validity of any modelling of the flux function far from this 
maximum as it is sometimes done in the literature. 
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